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ABSTRACT 



In this paper we consider the contributions of anomalous commutators to var- 
ious QCD sum rules. Using a combination of the BJL limit with the operator 
product expansion the results are presented in terms of the vacuum condensates 
of gauge invariant operators. It is demonstrated that the anomalous contributions 
are no negligible and reconcile various apparently contradictory calculations. 



1. Introduction 



The use of canonical commutators in the evaluation of current algebra re- 
lations has produced many results whose effects are directly measurable. Still in 
many cases the canonical evaluation of the commutators is ill defined, as clearly ex- 
emplified by Schwinger's calculation 1 of (0 | [j°(0, x), J l (0, y)] | 0) for a conserved 
current J^. If a fermionic current of the type vfj^^if) is replaced in the above expres- 
sion, and canonical commutation relations are used the above expression vanishes. 
In contrast, using general principles (such as Lorentz covariance and the absence 
of negative norm gauge-invariant states), the above vacuum expectation value is 
seen to be non-zero. 

This fact has been used repeatedly (although somewhat sporadically) in the 
calculations of anomalous (i.e., non-canonical) contributions to various commuta- 
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tors, especially in the context of anomalous theories . Similar effects have also 
been shown to modify the (canonically obtained) properties of the electroproduc- 
tion sum rules 3 . 

Faced with these problems in the canonical evaluation of commutation relations 
an alternative definition of the commutators was proposed by Bjorken and by 
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Johnson and Low . This definition preserves all the desirable features of the theory, 
is well defined and coincides with the canonical results whenever the latter are also 
well defined. 

The starting point of the Bjorken- Johnson-Low (BJL) definition of the com- 
mutator of two operators A and B is the time ordered product T(AB), presented 
as a function of the momentum transfer p. One then obtains the Laurent expan- 
sion of this operator in p° (the energy transfer). The term proportional to l/p° is 
identified as (the Fourier transform of) the equal time commutator Terms in this 
expansion containing positive powers of p° are associated with the covariantizing 
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of the time-ordered product and can be ignored. 
#1 It is of course possible for this quantity to be divergent. 
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The applications of this method have been largely restricted to perturbation 
theory (see however Refs. 6,7). On the other hand, many interesting applications 
of current algebra reside in the area where perturbation theory cannot be applied. 
In order to use the BJL definition in a wider range of situations we first note that 
the commutator is obtained by studying the relevant time-ordered products in the 
limit of large energy transfers, and, therefore, that an operator product expansion 8 
(OPE) is appropriate. The procedure which we follow is therefore to perform an 
OPE of the said time-ordered product, to then use renormalization group argu- 
ments to determine the high-energy behavior of the coefficient functions and thus 
extract the terms that contribute to the commutator. The result is then expressed 
in terms of the residues of the coefficient functions multiplied by the matrix el- 
ements of the local operators appearing in the OPE . In these calculations all 
symmetries of the theory are manifest, and so the resulting commutator will also 
respect them. A similar method was proposed by Crewther many years ago 6 but 
was not developed significantly. 

This paper is organized as follows. In the following section we describe the 
method in detail. Section 3 presents a comparison of the present method with 
some explicit perturbative calculations in 1+1 dimensions. In sections 4 and 5 we 
consider the anomalous commutator modifications to the current algebra approach 
to the U(1)a problem. Section 6 presents some explicit calculations pertaining the 
general arguments presented in sections 4 and 5. Parting comments are presented 
in section 7. The appendix contains the comparison of the present method with 
the results of perturbation theory for a 3+1 dimensional model. 



#2 This is reminiscent of the results obtained using sum rules . 
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2. Description of the method 



In this section we develop a useful technique which allows us to extract in- 
formation about non-canonical contributions to equal time commutators without 
going through the lengthy steps of loop calculations involving triangle, box and 
even pentagon diagrams. This technique, which does not rely on perturbation the- 
ory, is based on the Bjorken- Johnson-Low definition of the equal time commutators 
and on the operator product expansion (OPE). 

According to the BJL limit prescription, the definition of the commutator of 
two operators is obtained from the high energy behavior of Green's functions 



lim i P0 [ d n x e~ ip - x (a \TA{x/2)B{-x/2)\f3) 

po^oo J 



= I d»-WP-*<a|[A(0,x/2) , 5(0,-x/2)]|/3>, 



(2.1) 



where po is the time-like component of the four momentum . The time ordered 
product T is not Lorentz invariant and differs from the corresponding covariant 
Green's function by terms involving delta functions of xq and its derivatives (cor- 
responding to a polynomial in po in momentum space). If the left hand side of 
(2.1) is evaluated using covariant perturbation theory then all polynomials in po 
should be dropped. The covariant time-ordered product will be denoted by T* 

Since we are interested in the large momentum behaviour of (a\T*AB\(3), it 
is appropriate to express this object as a sum of local operators (OPE) where the 
coefficient functions summarize the po — > oo behaviour, 

/ d n x e~ ip - x (a\T*A(x/2)B(-x/2)\ (3) = J2 C AP) (<x\O r \P), (2.2) 

J r 

where the local operators O r are evaluated at x — 0. Taking the BJL limit of the 



#3 Double commutators can be denned using a straightforward generalization involving a dou- 
ble limit. 
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previous expression we find 




in— 1 



x e 



■ x (a|[A(0,x/2), £(0,-x/2)]|/3> 



(2.3) 



Elim [ipocr(p)] (a|0 r |/3), 



r 



where, as mentioned above, all terms in the c r growing as a power of po can be 
dropped. 

It is well known that the matching of dimensions of the operators T(AB) and 
O r in the OPE must take into account the anomalous dimensions of these objects. 
This can be avoided provided the operators considered are renormalization group 
invariant, such as the trace of the energy momentum tensor or the fermion mass 
terms. Note that even if A and B are renormalization group invariant, the time 
ordered product T(AB) need not have this property. In the most favorable cases 
the operators are renormalization group invariant and the canonical evaluation of 
the dimensions remains valid. 

Another characteristic of the method is that the results are evaluated in terms 
of a set of unknown constants, the residues of the the coefficient functions c r . For 
the applications which we consider this will not be a disadvantage: these constants 
multiply the matrix elements [a\O r \f3) which, in most cases cannot be evaluated 
to all orders in perturbation theory. Thus the final result will be given in terms of 
these "condensates" multiplied by the said constants. 
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3. Simple example 



As an application the previous remarks we consider a model containing fermions 
coupled to external non-Abelian gauge fields. We then choose A = J®, B — 
where J denotes the right or left-handed, gauge invariant fermionic current, and 
a, b denote color indices. Thus we consider 



q-ab 
1 [IV 



J d n x e~ ip - x T*J«{x/2)Jl 



-x/2). 



Similarly we define 



-tab 



lim ip T' 



ab 



(3-1) 



(3.2) 



(where the terms growing like a polynomial in po are dropped as discussed in the 
previous section). In writing the operator expansion of this object we have to 
pick terms that have the same dimensionality and that posses same symmetries; 
in particular we can restrict ourselves to gauge invariant operators. 

In n dimensions T ~ (mass ) n 2 , hence we can restrict ourselves to operators 
O r of dimension equal or greater than (n — 2) on the right hand side of the OPE. We 
will consider here the 1+1 dimensional case, leaving the 1+3 case to the appendix, 
as it does not bring up any new ideas or physics, and it is somewhat more involved. 

In 1+1 dimensions T has canonical dimension zero. Moreover, the Dirac ma- 
trices satisfy 757^ = e^ Vr y v which implies that we need only consider the vector 
currents, which we denote by J®. The current is given explicitly by 



/ x+e/2 



V exp 



\ 



A a (y)dy 



\c-e/2 



/ 



lb 



(x-0+h.c., (3.3) 



# 4 «i 



where T a denote the (anti-hermitian) group generators, " "h.c" denotes the Her- 
mitian conjugate, and = AaT a denotes the gauge field (the coupling constant 
is absorbed in the definition of A). 



#4 The conventions we use are tr(T a T b ) = -S ab , [T a ,T b ] = f abc T c , trT a {T b ,T c } = -d a 



be- 
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The OPE of T (defined now for the vector currents by replacing J — > J in 
(3.1)) is then given by 

T$ = co^l + ci ^(P)^ + c 2 % p {p)JP + • • • , (3.4) 

where 1 denotes the unit operator, J p is defined in (3.3), J p denotes the singlet 
vector current ~ iplpip, the ellipsis denote terms which will not contribute once 
the BJL limit is taken ; the c- number functions c r must have dimension — 1 

for r = 1,2 and for r = 0. Using the fact that T must be symmetric under 

#6 

a <-> b, [i <-> u, p <-> — p, the most general form of the coefficient functions is 



1 

p 2 



c W(p) =—dabc [a {g W Vv - gvpPn) + (3 {tppVv - tvpVp) + l^pvPp] 



^ 9 fabc 



| (#mp^ + 9u P Pp) + ^ (e w p„ + tvpVp) ( 3 5 ) 

, ,PfiPuPp 

+cgpvPp + a ^ 2 

c 2pt p (p) = p \ s ab [u {gppPv - g vp Pp) + v {e^v - e vp p p ) + we pu p p ] . 

Denoting by P the spatial component of the momentum the commutator for the 
vector currents, obtained by replacing J — > J in (3.1) and (3.2), is given by 

-C$=f abc [(a + c + d)4 + bJt}, 

% 

±C$ = - V^abP ~ d abc [(P + 7 ) J C + «^i C ] + \fabc K + aJt) 

~ ^ab [(v + w)J + uJi] , 

-Cfi = — cfafcjQ. 

I 

In the case where the theory has only right-handed couplings these relations im- 



#5 Note that the coefficient functions for the operators ipT c ip, ijiip, i/;T c j 5 ip and ^75^ will be 
of the form ~ (mass parameter)/}? 2 and so will not contribute in the BJL limit. 

#6 The constants £, . . . , w can be evaluated perturbatively. We will not need their explicit 
expressions. 

#7 In our conventions e 01 = +1. 
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ply d — and 



i[j a R (0,X/2),J b R (0,-X/2) 

This expression should be compared to those obtained in Ref. 10 which has the 
same form, except that Jr is replaced by Ar = Aq + A\. The discrepancy can 
be understood by following the procedure used in Ref. 10. What was done was 
to evaluate various matrix elements of the commutators and then to exhibit some 
local operators which have the same matrix elements. These operators are not 
unique, however. For example the matrix elements of Ar and Jr between the 
vacuum and the one gauge-boson state are proportional to each other in the zero 
momentum limit (the limit in the case of Jr is taken symmetrically, first averaging 
over the direction of the momentum and then letting the magnitude go to zero). It 
is easy to see that the results of the diagrammatic calculations are consistent with 
those presented in Ref.10 when Ar is replaced by —2hJr. It is in this sense that 
the above calculation in consistent with the explicit diagrammatic evaluation (up 
to the undetermined constants i] and a + b which we do not evaluate at this point.) 
We also point out that the above expressions have the expected form when taking 
the matrix elements of the commutators for states containing fermions. 

The above expressions of the anomalous commutators have the additional ad- 
vantage of being manifestly gauge covariant. The terms proportional to 8(X) are 
generated by the matrix elements of the canonical contribution to the commutator. 
The only irreducible non-canonical contribution is the Schwinger term oc 5^5' (X). 
We shall see in the appendix that similar results hold in 3+1 dimensions. 

The above results can also be used for calculating the Schwinger term and 
seagull for the commutators under consideration. Writing the expressions using a 



#8 Note that in 1+1 dimensions J^ M= o = Jrh=i = Jr 



= -iSabS'(X) - ^-f ahc J c R 8{X) (3.7) 
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time-like unit vector n we obtain for the commutator of two vector currents, 



J a R (x/2), J b R (-x/2)} 5(x ■ n) = Cf^ 2 \x) + S# a d a 5^ 2 \x) 



where Cf v = -f abc (\g^\ n ■ J c + \e^\ 
S$ a = 6 ab r t \e l „,\n p e'> a 



pa J c o) ; 



(3.8) 



From this it follows that the corresponding seagull vanishes. 

Thus the method is seen to work to lowest order in perturbation theory. The 
disadvantage is that the final result is expressed in terms of a few unknown con- 
stants which, if required, can only be obtained doing detailed calculations. More- 
over, for higher orders in perturbative calculations the anomalous dimensions of the 
various operators must be taken into account. We have seen that the apparently 
gauge variant results obtained in the literature can be re-interpreted as generated 
by the canonical terms in the commutator. 



4. Current algebra and the U(l) problem. 

In this section we will consider the effects of anomalous commutators in the 
study of the U(1)a problem. In this area the results obtained using instanton 

11 12 

calculations were criticized on the basis of certain inconsistencies which arise 
when the commutators involved are evaluated using canonical expressions. We will 
see that the relations derived in Ref. 12 are in general modified due to the anoma- 
lous terms in the commutators; this point is also made in Refs. 13 and 14 where 
it is noted that configurations carrying topological charge affect the pion decay 
constant Reference 12 also points out several apparent contradictions concerning 
the periodicity of the angle within the instanton and the canonical approaches. 
This problem was investigated in Refs. 13 and 14 and found to be rooted in a 
mis-application of the index theorem for which there are subtleties connected with 
spontaneous symmetry breaking. 
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These modifications are sufficient to explain the differences between the two 
approaches. 

To specify the notation we denote by J5 the gauge-invariant anomalous current 
which, in the presence of i massless flavors, satisfies 

dpJg = tv; u = ^F-F. (4.1) 

The charge associated with this current is denoted by Q5 = / d 3 x J5. 

In the effective lagrangian description the effects of an instanton (respectively, 
an anti-instanton) localized at x is described by a potential 14 U a (respectively, U*) 
with the identification (£ is the number of light quark flavors) 

£ u(x) ^81mU a (x), (4.2) 

where v is defined in (4.1). The potential U a is proportional to a quark-determinant 
operator involving all light flavors 11 ' 14 . 

The problem arises because U & has chirality 21 and so the right hand side of 
(4.2) has non-trivial commutator with the axial charge (as constructed in the effec- 
tive theory). In contrast, v apparently commutes with Q5, thus raising questions 
about the above identification. This contradiction can be solved by using the BJL 
definition of the commutator between v and J§. 

As a fist step we consider, for example, the vacuum correlator 

/ d A x e- ip - x (0\T*v(x/2)Jg(-x/2)\0) = J2 c r(p) (0 |O r | 0). (4.3) 

J r 

The lowest dimensional (non-trivial) operator that contributes to the right hand 
side is the trace of the energy- momentum tensor which we denote by 0. We 
expect the commutator to be a renormalization group invariant quantity; in this 
case the coefficient function associated with will have the form Cq(p) = CQp^/p 2 . 
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We hasten to point out that the OPE is valid only for large p and so one cannot 
interpret this form of Cq(p) as corresponding to a massless pseudoscalar excitation. 
We then obtain 

J rf 4 xe ip ' x <0 | [i/(0,x/2) , J 5 °(0, -x/2)] | 0> = % c@ (0 |Oj 0). (4.4) 

Since we expect both cq and (0 |0| 0) to be non-vanishing, it follows that the 
commutator of v with Q§ is non-trivial also within the context of QCD 

When quark masses are included the above equation is modified since more 
operators become available. Specifically, one can include on the right hand side of 
(4.3) a term containing the operator D, where 

D = 2j2™ f qfqf, (4.5) 
/ 

(rrif and gj denote the mass and field associated with the quark of flavor /). In 
this case the corresponding coefficient function in the OPE (4.3) tales the form 
Cj) — cdP^/p 2 , and (4.4) becomes 

J d 4 xe ip - x (0 | [i/(0, x/2) , j£(0, -x/2)] | 0) = i c® (0 jOj 0) + ic D (0 \D\ 0). (4.6) 
As we will see below, cq ; d do not in general vanish. Hence the commutator 

14 

between v and Q$ receives non-canonical contributions. Model calculations also 
show that the expression for the said commutator acquires a non-canonical piece 
proportional to (0 |@| 0). The right hand side of (4.6) should vanish for massless 
quarks in the p — > limit; this is verified within a specific model in section 6. 



#9 A straightforward perturbative calculation shows that, at least to one loop, c® ^ 0; see 
below. 
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5. Anomalous Ward identities 



In the previous section we remarked that the operator v can have non-zero 
commutator with the gauge invariant axial current Jg ; in particular [u, Q§] ^ 
0. These results are supported by a straightforward application of the effective 
lagrangian proposed in Ref. 14. It is of course possible for the constant cq (and cp 
is rrif 7^ 0) to vanish, but this would not be consistent with the effective lagrangian 
approach. We also point out that v will mix under renormalization with operators 
which have non-zero chirality. 

Should the above commutator be different form zero, the anomalous Ward 
identities will be modified. Consider then a gauge invariant operator O and define 

4° \p) = J d 4 x e-* x (0 | T J 5fl (x) O(0)| 0). (5.1) 



12 9 

The requirement that there be no light isosinglet pseudoscalars ' implies that 
p ■ will vanish as p — > 0. It follows that, by the definition of the T symbol, 



d 4 x (0 | T d ■ J 5 (x) O(0)| 0) + (0 \[Q 5 , 0(0)] | 0) 

d 4 x (0 | T A(x) O(0)| 0) + J d 4 x (0 | T iv{x) O(0) | 0) + (0 \[Q 5 , O(0)] | 0) 

(5.2) 



where 



A = 2i ^2 m f1f 75 <?/ (5.3) 
/=i 

and v is defined in (4.1); we have assumed that the anomaly equation, d • J5 = 
A + £1/, is an operator identity. 

Now, following Refs. 12 and 9, we consider (5.2) for the cases O = Iv and O = 
A; cancelling the correlator of A and v which appears in both these expressions 
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we obtain 

e 2 J d 4 x (o\Tu( x ) i/(o)| o> = (o|[g 5 , a(o) -eu(o)]\o)+ J d 4 x(o \ t a{x) a(o)|o). 

(5.4) 

The T-product on the right hand side equals the corresponding T* prod- 
uct. This is because the approach described in section 2 shows that there are 
no Schwinger terms in the equal-time commutator of A(x) and A(y); the corre- 
sponding seagull is therefore 
D = 2^2 mqq; we will write 



5 

sponding seagull is therefore zero . The commutator [Q5, , A] is proportional to 



i[Q 5 , A] =2(1 + 5) D, D = 2j2m f q f q f (5.5) 

/=1 

where 5 = if the commutator is evaluated canonically. Finally we have i[Qs, v\ = 
cqQ + cjjD, where 6 denotes the trace of the energy-momentum tensor. Thus we 
obtain, for the case of two light flavors (£ — 2), 

m u m d 

+ md) 2 
(5.6) 



2i [ d 4 x (0 I T u(x) i/(0)| 0) = (5-c D ) (0| J D|0)-c e (0|9|0)-4/ 7r m2-^ 



The fist two terms come from the anomalous commutators, while the last term is 
generated by the canonical commutator and the TAA contributions as evaluated 
in Ref. 9. 

The above calculation show that in general we can expect deviations form 
the canonical expression for the dependence of measurable quantities on the CP- 
violating parameter 9. It is of course possible for the non-canonical terms to van- 
ish, still explicit perturbative calculations and effect ive-lagrangian arguments favor 
c&,D, 8 7^ 0. Since the dependence on 9 disappears from all physical observables 
when one fermion is massless, we expect ce(0 |@| 0) — 5(0 \D \ 0) to vanish when 
any quark mass is zero. We can then write 

i J d*x (0 I T u(x) i/(0)| 0) = -2f 2 ml(l - A) ( ^"^ )2 , (5-7) 

for some constant A. The conditions under which A = 1 (or even if this is at all 
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possible) cannot be determined using general arguments. We will see below that 
low-energy models of the strong interactions predict 1 — A = 0(1) (see below) so 
that the estimates of physical quantities on the trong CP angle are altered only 
by a factor of 0(1) (except, of course, in the case A = 1). 

In the following section we present several calculations where the various coef- 
ficients and vacuum condensates are evaluated within explicit models. 

5.1. Explicit calculations 

In this section we present several computations. We have evaluated the OPE 
coefficients cq, cjj and 5 to lowest order in QCD; we also evaluate the condensates 
(0 |0| 0) and (0 \ D\ 0) in a chiral model of the strong interactions. 

Perturbative calculations 

We first consider the calculation of Ca- The OPE of the product T*Jj*v 
contains, to lowest order, three operators: D defined in (4.5), O^, the energy- 
momentum tensor, and O, its trace. For the calculation at hand we evaluate the 
matrix element between the vacuum and a two gluon state, the relevant graphs are 
presented in figure 1. 



p-q/2 




i k + crossed 



Fig. 1 Lowest order diagrams contributing to the OPE coefficient of 
6 in the operator T*J^v. 

In order to obtain the OPE coefficients in the p — > oo limit we need only 
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consider the diagrams contracted with the momentum carried by the axial current, 
i.e., we multiply each graph by (p + q/2) fl (the momentum of the axial current) 
and contract the index \i. A simple calculation using 



= ^fc) 2 + 2-(l + 7m)A (5.8) 

where (3 denotes the beta function for a = g 2 /An and 7 m the mass anomalous 
dimension. Assuming the presence of I fermions we obtain 

c e = ^--; P(a) = - U ~ 2i/3 a 2 + 0(a 3 ). (5.9) 



We now consider cjj which can be extracted from the matrix element of TvJ§ 
between the vacuum and a two fermion state. It is easy to see, however, that the 
graphs for this matrix element are all 0(a 2 ); in contrast cq = 0(a). The OPE for 
TuJ^ contains both 6, defined in (5.8), and D. It follows that the 0(a) term in 
cq should be cancelled by a similar term in cd, hence 

c D = \cq (5.10) 

Finally we calculate 5 by evaluating the matrix element of T*A between the 
vacuum and a two fermion state. We skip for brevity the description of the graphs; 
the final result is 

to first order in a. 

these calculations show that, at lest to first order in perturbation theory, all 
the anomalous coefficients are finite and non- vanishing as claimed previously. 
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5.2. Model Calculations 

In order to obtain estimates of the condensates (0 |@| 0) and (0 \ D\ 0) we con- 
sider a chiral model which obeys the same symmetries as QCD. In order to gener- 
ate Green's functions involving v we modify the QCD lagrangian by adding a term 
9u/2, with 9 an external source (for details see Ref. 15). The Lagrangian takes 
the form 15 

C = -V + Vitrd^d^U + (V 2 trMU + H.c) + V 3 d(f) -d9 + V 4 (09) 2 (5.12) 

where Vi = l^(0o + 9), re& l? an d Vi(a) = V*(—a). The meson field, denoted 
by U, belongs to the unitary field U(£); we will write U = exp(i<j)o/£) E with 
£ G SU (£) . For £ = 3 £ describer the pseudoscalar meson octet and 0o describes 
the pseudoscalar isosinglet (the r/). 

The field S describes the usual pseudoscalar meson multiplet (under SU(£) 
flavor); 0o describes the pseudoscalar singlet {i.e., the r\ for t = 2, and the rj' for 
£ = 3). 

This model is an accurate representation of QCD at low momentum transfers, 
so we will not use it in obtaining the BJL definition of the commutators (which 
involve the p° — > oo limit). We can, however, use this model to evaluate the 
condensates (0 |0| 0) and (0 \D\ 0) and the low momentum limit of T*uu. We will, 
for simplicity, work in the SU(£) symmetric limit where M = ml. 

To lowest order in a momentum expansion the correlator (0 \T\ 0)*uu can be 

15 

obtained by replacing U in C by the solution to the classical equations of motion . 
A simple calculation shows that the lagrangian then takes the form 

£ = \9K9 + const.; K = -2 [y 4 (0) + invW^O) - V 3 (0)} - —V 2 (0) 

(5.13) 

which shows that the vacuum correlator (0 \T*vv\ 0) is proportional to m in the 
limit of zero momentum transfer. 

The various condensates can also be evaluated within this model. From Ref. 
15 we get 

(0 \D\ 0) = -Wmlfl- (0 |Oj 0) = 4Vb(0) - Mm 2 jl (5.14) 

where denotes the (degenerate) meson mass and f % the corresponding decay 
constant. Perturbative calculations suggest that cq remains non-zero as m — > 0, 
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hence consistency of the OPE with the above expressions for T*uu, (0 \D\ 0) and 
(O|0|O) requires Vq(0) to vanish as m — > 0. It follows that both condensates 
(0 \D\ 0), (0 |6| 0) vanish in this limit. This result justifies the claims made at the 
end of the previous section concerning the behaviour of the condensates in the zero 
mass limit. Similar results are obtained using the (closely related) model of Ref. 14 

Within this model i J d 4 x(0 \T*vv\ 0) = —^(m^f^) 2 which corresponds to A = 
in (5.7). Thus the possibility of having A = 1 and a dynamical cancellation of 
the dependence on the strong CP angle is not realized, at least within this model. 

6. Conclusions 

In this paper we considered the BJL definition for the commutators and applied 
it in conjunction with the operator product expansion. The method can be applied 
both in the perturbative and non-perturbative regimes. As an application of the 
first case we considered the anomalous commutator between chiral currents in 1+1 
and 3+1 dimensions. We showed that in this case the results in the literature can 
be re-interpreted to yield a gauge invariant expression for the commutators. The 
method here proposed is consistent with these results. 

In the non-perturbative regime we considered the current algebra relations be- 
tween the instanton number density v and the gauge-invariant anomalous axial 
charge. We showed that, in general, this commutator is no n- vanishing, in accor- 
dance with the results obtained using instanton calculations. We also noted that 
this conclusion is based on the no-trivial chiral transformation properties of the 
instanton density and this leads to some modifications of the expressions resulting 
from the anomalous Ward identities. 

The method requires some knowledge about behaviour of the coefficient func- 
tions (which appear in the OPE) at large momentum transfers p. In asymptotically 
free theories this is available via the renormalization group. The final results are 
expressed in terms of the residues of the coefficient functions (i.e., the constant 
multiplying the term behaving as 1/po) and of the matrix elements of various local 
operators (the "condensates"). These quantities can be evaluated explicitly within 
perturbation theory; in the non-perturbative regime the condensates cannot be 
evaluated explicitly but can be used to parametrize the results. 

Whereas the OPE coefficients can be evaluated perturbatively to any desired 
order of accuracy, the condensates are not calculable in this manner; for these 
quantities effective models must be considered. Unfortunately the effective theories 
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are valid only at small momenta and this implies that the p — > oo limit of the OPE 
coefficient functions cannot be accurately evaluated using these theories. Explicit 
calculations verify several claims made on general grounds: there are non-trivial 
non-canonical contributions to the commutators. These contributions can be used 
to reconcile the operator and instanton approaches to the U(1)a problem. 
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APPENDIX 

In this appendix we consider the more complicated case of the OPE of the 
two current correlator in 3+1 dimensions. In order to keep the discussion at a 
manageable level we will consider the case 

T^( P ) = J d*x e-^T*J^(x/2)4(-x/2). (A.l) 

Now T has mass dimension = 2, leading us to a more complicated OPE. The 
relevant terms are 7 ^ 10 



T ab (v) = — 



VOL 



where 



+ («? "DA + ufV v F^ + uf c V a F^) p a 
+ < bc J c ^ Pv + v#*J$ vPlt + uf c J C YW 

+ Kb {vih^Pv + v lhvPii + V< i J p p G e puVpa ) 



uf c = u^f abc + v!fd ahc 



(note that V a F afl = due to the Bianchi identities). 



(A.2) 



(A.3) 



#10 We use the following conventions: F° v = \t^ p „FP\ with F c pa = d p A% - 8 a A c p + f abc A a p A b a . 
The covariant derivative is given by V p F c pa = d^F^ + f dec A d p F^. 
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Applying the BJL limit to the previous expression, we obtain 



lim Po72?(p) 



d 3 xe ip -* 



J« 5 (0,x/2),4(0,-x/2) 



u 



abc 



fj,0 



9viKo + 9riK ) + uf c ( g^Fto + g^F t 



P' 



(A.4) 



The equal time commutator for the ji — v — case is given by 



J a 5 (0,x/2),J b (0,-x/2) 

= -i (uf c + uf c ^ B c • V<J(x) 



(A.5) 



+ 



(uf c + c ) J 5 C + <U (ui + v 2 ) J50] 5(x), 



where B c denotes the chromo-magnetic field, B l = F^. The commutator for the 
space component of the vector current and the time component of the axial current 
is 



J a 5 (0,x/2),J b (0,-x/2 

= iuf C E C x V5(x 



(A.6) 



+ 



{uf c - uf c ) V B C + uf C 3 C 5 + 8 ab v 2 J 5 ] 5(x), 



where E c denotes the chromo-electric field, E % = F^. 

These results are, as in the 1+1 dimensional case, manifestly gauge covariant. 
We have verified that they are consistent with the explicit loop calculations pre- 
sented in ; for example the two-gauge-boson matrix elements of id^efeadF^ A c k 
and M 2 f abc J c R ^ =Q coincide. 

These results can be used to calculate the corresponding Schwinger terms and 
covariantizing seagulls. Following the procedure described in Ref. 5 we obtain 
[J% 5 (x/2),4]5(n ■ x) = Cf v 5^\x) + S^.^d^ix) where the Schwinger term 
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equals (the computation is straightforward and only the results will be presented) 

S;l ;a = - [uf c (g va Ffo + g vp F; a ) + uf c (g^FSp + g^F° a ) ] nP ( A.7) 

and the corresponding seagull is 

= - {< bc n v F; a + uf c n^ a ) n a (A.8) 

We also remark that (again following the procedure described in Ref. 5) when both 
currents are conserved the requirement that the T* product be Lorentz covariant 
implies uf c = uf c , uf c + uf c + uf c = and uf c = uf c = v x = v 2 = 0. Finally 
we note that the conditions under which the Schwinger terms cancel against the 
seagull contributions to the Ward identities is simply uf c + uf c = 0. 
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